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Abstract 

In this paper we introduce a simple model for a financial market characterized by a sin- 
gle stock or good and an interplay between two different traders populations, chartists and 
fundamentalists, which determine the price dynamic of the stock. The model has been in- 
spired by the microscopic Lux-Marchesi model |18] . The introduction of kinetic equations 
permits to study the asymptotic behavior of the investments and the price distributions and 
to characterize the regimes of lognormal behavior and the formation of power law tails. 
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1 Introduction 

Most speculative markets at national and international level share a number of stylized facts, like 
volatility clustering and fat tails of returns, for which a satisfactory explanation is still lacking in 
standard theories of financial markets [2S] . Such stylized facts are now almost universally accepted 
among economists and physicists and it is now clear that financial markets dynamics give rise to 
some kind of universal scaling laws. 

Showing similarities with scaling laws for other systems with many interacting particles, a 
description of financial markets as multi-agent interacting systems appeared to be a natural con- 
sequence [TBJ [TS1 [221 [551 [33J . This topic was pursued by quite a number of contributions appearing 
in both the physics and economics literature in recent years [H [31 HI 121 [HI (HI H51 1221 12S1 133J ■ This 
new research field borrows several methods and tools from classical statistical mechanics, where 
emerging complex behavior arises from relatively simple rules due to the interaction of a large 
number of components |24j . 

Starting from the microscopic dynamics, kinetic models can be derived with the tools of classical 
kinetic theory of fluids [H El El [SI H31 El 1201 [231 [2H ED] - In contrast with microscopic dynamics, 
where behavior often can be studied only empirically through computer simulations, kinetic models 
based on PDEs allow us to derive analytically general information on the model and its asymptotic 
behavior. 

In this paper we introduce a simple Boltzmann-like model for a speculative market characterized 
by a single stock and a socio-economical interplay between two different types of traders, chartists 
and fundamentalists. The model is strictly related to the microscopic Lux-Marchesi model |18) 
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and to kinetic models of opinion formation recently introduced in [32 . In addition, we take into 
account some psychological and behavioral components of the agents, like the way they interact 
each other and perceive the risk, which may produce non rational behaviors. This is done by means 
of a suitable "value function" in agreement with the Prospect Theory by Kahneman and Tversky 
[121 113) . As we will show people systematically overreacting produces substantial instabilities in 
the stock market. 

In an earlier paper [8] a similar approach has been used considering a single population of 
investors interacting in the stock market on the basis of the microscopic Levy-Levy-Solomon model 
[T5] . The emergence of a lognormal behavior for the wealth distribution of the agents has been 
shown. Though the theoretical set-up of the analysis is close in certain respects to that of [5] , the 
structure of the model is rather different. Namely, the description of individual behavior follows an 
opinion formation dynamic strictly connected with the price trend. In this way, the heterogeneity 
among agents as well as their social interactions will be taken into account which both are key 
elements affecting the outcome of the overall market dynamics. 

Following the analysis developed in [32] , we shall prove that the Boltzmann model converges 
in a suitable asymptotic limit towards convection-diffusion equations of Fokker-Planck type. Other 
Fokker-Planck equations were obtained using different approaches in [I] [3TJ [31] . This permits to 
study the asymptotic behavior of the investments and the price distributions and to characterize 
the regimes of lognormal behavior and the ones with power law tails. The main finding of the 
present paper is that the presence of heterogeneous strategies, both fundamentalists and chartists, 
is essential to achieve basic stylized fact like the presence of fat tails. 

The rest of the paper is organized as follows. In Section 2 we introduce the Boltzmann kinetic 
model for the interacting chartists and the price evolution. Details of the strategy exchange 
between chartists and fundamentalists are also presented here. A characterization of the admissible 
equilibrium states of the resulting system is then reported. Next, in Section 3, with the aim 
to study the asymptotic behavior of the chartists and price distributions, we introduce simpler 
Fokker-Planck approximations of the Boltzmann system and give explicit expressions of the long 
time behavior. The mathematical details of the derivation of such Fokker-Planck models are 
reported in separate appendices at the end of the manuscript. Numerical results which confirm 
the theoretical analysis are given in Section 4 and some concluding remarks are discussed in the 
last section. 

2 A kinetic model for multiple agents interactions 

We describe a simple financial market characterized by a single stock or good and an interplay 
between two different traders populations, chartists and fundamentalists, which determine the price 
dynamic of such stock (good). The aim is to introduce a kinetic description both for the behavior 
of the microscopic agents and for the price, and then to exploit the tools given by kinetic theory 
to get more insight about the way the microscopic dynamic of each trading agent can influence 
the evolution of the price, and be responsible of the appearance of 'stylized' fact like 'fat tails' and 
'lognormal' behavior. 

2.1 Kinetic setting 

Similarly to Lux and Marchesi model |18j , the starting point is a population of two different kind of 
traders, chartists and fundamentalists. Chartists are characterized by their number density pc and 
the investment propensity (or opinion index) y of a single agent whereas fundamentalists appear 
only through their number density pp- The value p = pf + Pc is invariant in time so that the 
total number of agents remains constant. In the sequel we will assume for simplicity p = 1. 
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Dynamic of investment propensity among chartists. Let us define f(y,t), y E [—1,1], 
the distribution function of chartists with investment propensity y at time t. Positive values of y 
represent buyers, negative values characterize sellers and close to y = we have undecided agents. 
Clearly 



PC® = jf* 



f(y,t)dy. 



Moreover we define the mean investment propensity 

Y(t) = ^— f 1 f(y,t)ydy. (1) 

pc{t) J-i 

For a given price S(t) and price derivative S(t) = dS(t)/dt the microscopic dynamic of the invest- 
ment propensity of chartists is characterized by the following binary interactions (y,y*) —> (y',y*) 
with 

?/ = (l-a 1 H(y)-a 2 )y + a 1 H(y)y* + a 2 <S> [ J + D(y)r), 



S(t) 



y'* = (l-a 1 H(y*)-a 2 )y* + a 1 H(y*)y + a 2 i f> (f|yj + D {v*)V* 



(2) 



Here a,\ E [0, 1] and a 2 E [0, 1], with a\ + a 2 < 1, measure the importance the individuals place 
on others opinions and actual price trend in forming expectations about future price changes. The 
random variables rj and 77* are assumed distributed accordingly to Q(rj) with zero mean and variance 
<j and measure individual deviations from the average behavior. The function H (y) E [0, 1] is 
taken symmetric on the interval /, and characterize the herding behavior, whereas D(y) defines 
the diffusive behavior, and will be also taken symmetric on /. Simple examples of herding function 
and diffusion function are given by 

H(y)=a + b(l-\y\), D{y) = (l-y 2 )\ 

with 0<o + 6<l, o>0, 6>0, 7>0 (see Figure f27Tj) . Other choices are of course possible, 
note that in order to preserve the bounds for y it is essential that D{y) vanishes in y = ±1. Both 
functions take into account that extremal positions suffer less herding and fluctuations. For 6 = 0, 
H(y) is constant and no herding effect is present and the mean investment propensity is preserved 
when the market influence is neglected (a 2 = 0) as in classical opinion models a model ( see [52] 
at the reference therein). 
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Figure 1: Typical examples of herding function H(y) (left) and diffusion function D(y) (right). 

A remarkable feature of the above relations is the presence of the normalized value function 
$>(S(t)/ S(t)) in [—1, 1] in the sense of Kahneman and Tversky [12l [13] that models the reaction 
of individuals towards potential gain and losses in the market [12] . This permits to introduce 
behavioral aspects in the market dynamic and to take into account the influence of psychology on 
the behavior of financial practitioners. 

The value function is defined on deviations from a reference point, which is usually assumed 
equal to zero (but it can be considered also positive or negative) , and is normally concave for gains 
(implying risk aversion), commonly convex for losses (risk seeking) and is generally steeper for 
losses than for gains (loss aversion) (see Figure 12. 1|) . 



Figure 2: An example of value function $>(S(t)/S(t)). 

Let us ignore for the moment the price evolution. The above binary interaction gives the 
following kinetic equation for the time evolution of chartists 

% = Q(f,f), (3) 

where for any test function ip the interaction operator Q can be conveniently written in weak form 
as 

/ Qv>{y)dy= / B(y,y*,r),r)*)f(y)f(y*)(tp(y') -<p(y))dr]dr)*dy*dy (4) 

J-i J[-i,h 2 Jb 2 
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with the transition rate has given by 

being %(•) the indicator function. Note that the mass density of chartists pc(t) is an invariant for 
the interaction, (ip = 1). 

It is worth to observe that for a given Dc (y) a suitable choice of the support of the random variable 
r], avoids the dependence of the collisional kernel B(y,y*) on the variables y,y*. 
As an example, if we take D(y) = 1 — y 2 we have 

?/ = (l-a 1 H(y)-a 2 )y + a 1 H(y)y,+a 2 ^(^^ ] j +{l-y 2 )r 1 

< (1 - a 1 H(y) - a 2 )y + a\H(y) + a 2 + (1 - y 2 )v- 
Then to have y' < 1 for any y e [— 1, 1], wc have to chose j) such that 

{l-y 2 )ri< (1 - ai - a 2 )(l - y) 

which gives 

Analogously we can ensure y' > — 1 , thus it is enough to take 

e [-^(1 - ai - a 2 ), ^(1 - ai - a 2 )]- 
For this reason, in the rest of the paper, we will consider only kernel of "maxwellian type" 



B(y 7 y*,v,V*) = ©(»?)©(»?*). 



Strategy exchange chartists-fundamentalists. In addition to the change of investment propen- 
sity due to a balance between herding behavior and the price followers nature of chartists, the model 
includes the possibility that an agent changes its strategy from chartist to fundamentalists and 
viceversa. 

Agents meet individual from the other group, compare excess profits from both strategies and 
with a probability depending on the pay-off differential switch to the more successful strategy. 
When a chartist and a fundamentalist meet they characterize the success of a given strategy 
trough the profits earned by comparing 

Here ip(y) e [—1,1] has the same sign of y and takes into account the change of sign in the 
profits accordingly to the actual behavior of the agent in the market which rely on his investment 
propensity y. The simplest choice is ip(y) = sgn(y). 

The value D is the nominal dividend and r the average real return of the market, such that 
r = D/Sf, i-e. evaluated at its fundamental value Sf in a state of stable price S = the asset 
yield the same returns of other investments, or equivalently Xc = Xf = 0. The discount factor 
k < 1 is justified by the observation that Xp is an expected gain realized only after reversal to the 
fundamental value. Finally \i > measures the frequency of the exchange rates. 
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A chartist characterized by an investment propensity y and a fundamentalist meet each other, 
and after comparing their strategies, they exchange strategies with a rate given by a suitable 
monotone function B F c{-) > 0. More precisely a chartist switch to fundamentalist with a rate 
Bpc(Xp — Xc) and a fundamentalist switch to chartist at a rate Bpc{Xc — X F ). 

For chartists we define the following linear strategy exchange operator 

QfcU) = n PF {t)f{y){B F c{Xc - X F ) - B FC {X F - X c )) 

where p > measures the frequency of the exchange rates. 

Taking into account such strategy exchanges we have the chartists-fundamentalists model 

( df 

-£ = Q(f, /) + Wp(t)f(y)(B FC (X c - X F ) - B FC (X F - X c )) 
= pp F (t) J ^ f(y)(B FC (Xp - Xc) - Bpc(X c - X F )) dy. 
It is immediate to verify that the total number density pc + Pf is conserved in time. 



Price evolution. Finally we introduce the probability density V(s,t) of a given price s at time 
t. The effective market price S(t) is defined as the mean value 

/>oo 

S{t) = / V(s,t)sds. (7) 

Following Lux and Marchesi |18j the microscopic dynamic of the price is given by 

s' = s + l3(p c t c Y(t)s + PFl (S F - a)) + lis (8) 

where the parameters /3, represent the price speed evaluation, ry is a random variable with zero 
mean and variance £ 2 , distributed accordingly to 4 r (^). In the above relation chartists either buy 
or sell the same number to of units and 7 is the reaction strength of fundamentalists to deviations 
from the fundamental value. 

Thus the chartists-fundamentalists system of equations ^ is complemented with the equation 
for the price distribution 

dV 

where the operator L, is linear, and in weak form it reads 

L(V)(sMs)ds= b(8,rt)V(s)(<p{i/)-<p(s))drida (10) 



JO JM 

with the transition rate b(s,rj) = 4 , (fy)x( s ' > 0). 

As before, a suitable choice of the domain for the support of variable 77 ensures a' > 0. Assuming 

t] G [-1 + (3(pcT c + pfi), 1 - P(pcT c + Pfj)}, P{pcT c + pFl) < 1- 
permits to express the transition rate in the simpler form 

b{s,f]) = 9(rj). 

Note that the expected value for the stock price satisfies the same differential equation as 
in[12j[T8] 

dS(t) 



dt 



p Pc t c Y{t)S{t) + p PFl {S F - S{t)). (11) 



G 



Booms, crashes and macroscopic stationary states. In order to study the macroscopic 
steady states and relate them to the value function $ let us start by observing that the equilibrium 
states for the price satisfy 

Pc t c YS + Pf j{Sf -S) = 
and thus fall in one of the following categories 

(i) PF £ 0, S = PF P /J p S J cY , PFlS F - pct c Y > 0. 

(ii) pf = 0, Y = 0, S arbitrary, 

(iii) pf = 0, S = 0, Y arbitrary. 

At equilibrium we require pp, pc an d Y to be constants. In order for the number densities to be 
constants we require Qfc — 0. For pp ^ and pc ^= 0, thanks to monotonicity of Bpc, we have 
Xc — Xp or equivalently S — Sf. Note that Qfc vanishes also when pp = or pc = 0. These 
considerations reduce the set of possible equilibrium configurations to 

(i) pf^O, S = S F , Y = 0, 

(ii) pF = 0, Y = 0, S arbitrary, 

(iii) pp = 0, S = 0, Y arbitrary. 

Finally we consider the requirements for Y to be constant. In the case Qfc — the first moment 
equation reads 

J t Y (t)= - a 1 J^H(y)yf(y)dy-a 2 p c Y(t) 

's(t) 



+ 

which gives the steady state condition 



ai Y(t) J^H(y)f(y)dy + a2 p c $(^ 



»i J ^H(y)yf(y)dy-a 2 p c Y + a 1 Y J ^H(y)f(y)dy + a 2 p c $ = 0. 



This gives a constraint for the value function $, precisely 



a 2 p c <S> (jrl j = ai J H(y)yf(y)dy + a 2 p c Y - u x Y J H(y)f(y)dy 



which in the simple case of H constant reduces to 



Now using the fact that 



S{t) v/s {Sf - S(t)) 
we can state 
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Proposition 1 The system of equations (0|) in the case of H constant admits the following possible 
equilibrium configurations 



(i) 




s = 


Sf, 


Y = 0, 


$(0) = 0, 


(ii) 


PF = 0, 


Y = 


o, 


$(0) = 0, 


5* arbitrary, 


(Hi) 


PF = 0, 


Y = 


Y*, 


with = 


$(Pt c Y*), S = 



Note that if the reference point for the value function $(0) 7^ configuration (i) and (ii) are not 
possible for a constant H . This is in good agreement with the fact that an emotional perception 
of the market from the chartists acts as a source of instability for the market itself. In contrast 
configuration (Hi), corresponding to a market crash, can be achieved also for $(0) 7^ 0. The 
existence of a unique fixed point Y r has to be guaranteed by the choice of j3 and tc- Of course if 
the reference point is set to zero, $(0) = 0, we have Y* = 0. It is easy to verify that these possible 
equilibrium configurations include the ones in the original Lux-Marchesi model j!7j . 

In addition to the above equilibrium configurations the model admits several other possible 
asymptotic behavior in the form of booms and cycles. Some of the fundamental features of the 
model are summarized in the following. 

Remark 1 

• Chartists alone (pF = 0, pc = lj influence the price through their mean propensity to invest 
Y(t) and at the same time the price trend influences their mean propensity through the value 
function <$>(S(t)/ S(t)), since S(t)/S(t) — j3Y(t)tc- Thus, except for the particular shape 
of the value function, if the mean propensity is initially (sufficiently) positive then it will 
continue to grow together with the price and the opposite occurs if it is initially (sufficiently) 
negative. 

The market goes towards a boom (exponential grow of the price) or a crash (exponential decay 
of the price) with 

5(0)e^ fc < S(t) < S(O)e 0tc , 

and agents tend to concentrate in y = 1 and y = — 1 respectively depending on the choices of 
H and $ . This is in good agreement with the price followers nature of chartists. 

• Fundamentalists alone (pp = {,Pc = 0J influence the price through their expectation of the 
fundamental price. So their effect is to drive the price towards the fundamental price. For a 
constant fundamental price Sf the equilibrium state reached is characterized by S = Sf and 
the trend is exponential. 

• The presence of fundamentalists acts in contrast to the chartists pressure towards market 
booms or crashes. If their number is large enough they are capable to drive the price towards 
the fundamental value otherwise the chartists dynamic may dominate. In addition to booms 
and crashes, we have now the possibility of price cycles/oscillations around the fundamental 
value. 

3 Fokker-Planck approximations and asymptotic behavior 

Now we consider what happens at the kinetic scale. Due to the extreme difficulty to get detailed 
information on the asymptotic behavior of the kinetic coupled system, we will recover for both 
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distribution functions /, and V, simplified Fokkcr-Planck models which preserve the main features 
of the original kinetic model. To keep notations simple, since we are mostly interested in the study 
of the equilibrium states we ignore the presence of the terms describing the change of strategy. 
However they can be easily included in the scaling described below. 
For this purpose we introduce a time scaling parameter £ and define 



t = &, f(y,r) = f(y,t), V(s,r) = V(s,t). 
To preserve the chartists dynamic in the limit, we must require that 

o 

012 



lim — - 



Oil, 



lim 



£*2, 



lim — 



(12) 



(13) 



where A is a positive constant. 

Similarly for the price dynamic, we assume 



lim 4 

Performing similar computations as in 
following Fokker-Planck system 



lim — 



P, lim„ = v. (14) 
(see Appendix A and B for details) we recover the 



r 9[ 

dr 



d_ r 

dy 



—V + — 

dr ds 



PC (diH(y)(Y -y) + d 2 (# - y)) / 
P (pcYt c s + PFl (S F - a)) V 



Xp c d 2 
2 dy 



(15a) 
(15b) 



where we used the shorthand $ for <!> [S(t)/S(t) 

For notation simplicity in the sequel we will omit the tildes in the variables /, V, Y and S. 

If we now take D(y) = 1 - y 2 , and H(y) = 1 we can compute explicitly the equilibrium state 
for chartists with a constant mean investment propensity Y = Y* as 



_2+y.iS±£&l/, _,x_ 2 _ n (£i+£2) ( {1-Y t y)(di +d 2 y 



f°°(y) = C (l + y)- 2+y -^-(l-y)-'- r *S^e W [-i A1 _ y2) ) ( 16 ) 

where Co = Cq(Y*, A/ (a"i + 02)) is such that the mass of /°° is equal to pc- Other choices of the 
diffusion function originate different steady states (sec 32 ). 

Observe that, in the case Y* ^ 0, the distribution is not symmetric and in the chartist population 
a predominant behavior arise. Otherwise when the reference point of the value function is set to 
zero we have a symmetric distribution with two peaks and mean value zero, and the macroscopic 
state of indecision is given, microscopically, by a polarization of the chartist population among two 
opposite kind of behaviors (see Figure |3]). 

In order to study the asymptotic behavior for the price we must distinguish between the case 
Pf 7^ and pf = 0. 

Let us consider first the situation in which pf = (or equivalently pc = 1). For this purpose, 
we introduce the scaling 

V(s,t) = -v(x,t), x = 1 °g(s)- 
s 

It is straightforward to show that v(xi T ) satisfies the following linear convection diffusion equation 



1 



v - I d . v d 

2 " p c \ % u(x ' r) + 2 ap u(x ' r) ' 
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Figure 3: Equilibrium distribution function of the chartist investment propensity for different 
values of Y* = 0, 0.2, —0.2 (left) and corresponding behavior of the price S (right). Exact solutions 
with p = 1,0 = 0.1, t c = 1, A/(di + a 2 ) = 1 and f(y,0) = f°°(y). 



which admits the self-similar solution [8] 



, , 1 ( (x + log( V / ^)/^(r))-log(5(r))) 2 



(21og(£(T)/S(T) 2 )7r)5 ^ 2\og{E{r)/S\r) 
with 

E(t) = / V(s,T)s 2 ds. 
Jo 

Then reverting to the original variables it gives the lognormal behavior 



v( s 1 ( (l0g(s^MT)/S(T) 2 ? \ 

[S,T> S (2\og(E(r)/S(ry)n)i P { 2\o g (E(r) / S{rf) )' ( '> 
where E(r) satisfies the differential equation 

— = (2f3Yt c + v)E{T). 

Thus for a steady state characterized by (ii) in Proposition Q] we have S(t) — So, Y — and 
E(r)=e^E . 

Besides the above equilibrium state, equation (TT71) characterizes also the self-similar behavior 
of the price distribution in the case of booms and crashes, when the price S(t) grows arbitrary or 
decays to zero. In particular in the limit S(t) — > 0, point (Hi) in Proposition [TJ the distribution 
function V(s,t) concentrates near zero. 

Finally we consider the microscopic behavior of the model where both pc ^ and pp ^ 0. 

Recall now the Fokker-Planck equation for the price (|15bl) and consider the stationary case (i) 
in Proposition [T] The Fokker-Planck equation in such case reads 



dr ds 



v d 



P PFl (S F -s)V =-—(s 2 V). (18) 
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In this case the steady state can be computed as [TJ [7] and yelds 

1 ((i-l)Sir 

y-( s ) = C 1 ( M )^e r^, (19) 

where p — l + 2f3pF"f jv and Ci(/x) = ((/U— \)Sf') >1 j^ip) with T(-) being the usual Gamma function. 
Therefore the stationary state is described by a Gamma-like distribution with Pareto power law 
tails. 

Remark 2 

• The 'presence of fundamentalists is then essential in order to obtain fat tails in the price 
distribution. Their presence force the price to approach the mean value Sp in a way similar 
to the redistribution of wealth in the models proposed in fJl 0/. This feature seems to be 
essential for the development of power law behaviors. The stationary state for the price HSfy 
has in fact the same structure of the stationary states for the wealth in |IJ [^/. 

• In our description we have considered a constant value for the fundamental price. Such 
an assumption might seem quite unrealistic since, according to the economic literature, the 
fundamental price is usually treated like a temporal series with a stationary lognormal dis- 
tribution. This reflect the facts that the returns in logarithmic form are gaussian distributed 
with zero mean and a fixed variance, i.e big jumps between two successive realizations are 
rarely verified. Note however that introducing a given time dependent distribution function 
VsF(q,t) for the fundamental price such that 

r+oo 

S F (t) = / V S F(q,t)qdq, 
Jo 

and considering the following dynamic in the price evolution 

s' = s + (3(p c t c Y(t)s + p Fl {q - s)) + rjs, 
where s and q are random variables distributed as V(s,t) and VsF(q,t), we recover 

/>oo />oo />oo /> 

/ L(V)(s)<p{s)ds= / b(s,7])VsF(q)V(s)(ip{ S ')-<p(s))dr]dsdq 
Jo Jo Jo Jr 

which is the analogous of ilO\) and yields the same Fokker-Planck equation il5b\) for the 
asymptotic behavior of the model. We omit the details. 

4 Numerical examples 

In this section we considered different numerical simulations of the kinetic system. A Monte Carlo 
method analogous to the one used in classical rarefied gas dynamic has been used for the simulations 
[2]. In order to compute the kinetic behavior of the price, we use a set of N s = 50000 samples 
which can be though as possible realizations of the random variable s denoting the price. Since at 
the initial time the stock price So is supposed to be known, all samples are initialized at the same 
value initially. We represent the initial chartists distribution with a set of N c — pc(0)N sample 
agents with iV = 50000. These do not represent real agents but simply statistical realization of 
the random variable y. Such choices of N s and N permits to obtain results with a moderate effect 
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Figure 4: Equilibrium distribution function of the chartist investment propensity, with $(0) = 
(left) and log-normal distribution for the price (right) at t = 1500. The continuous line is the 
solution of the corresponding Fokker-Planck equation. 

of fluctuations without averaging. 

In all our computations we take the value function 



where x € [— L,L], Rq is the reference point and < I < r < 1. For example we choose r = 1/2 
and 1= 1/4. 



In the first test we consider the case with pp = i.e only chartists are present in the model. We 
computed the equilibrium distribution for $(0) = of the investment propensity. We take /3 — 0.1, 
tc — 1, a constant herding function H(y) = 1 and the coefficients oi\ = = 0.01. The initial 
data for the chartists is perfectly symmetric with Y = 0, so the price remains constant S = Sq 
with So = 10. A particular care is required in the simulation to keep Y = since the equilibrium 
point is unstable and as soon as Y ^ the results deviate towards a market boom or crash. 

After T = 1500 iteration the solution for the investment propensity has reached a stationary 
state and is plotted together with the solution of the Fokker-Planck limit in Figure |4] In the 
same figure we report also the computed solution for the price distribution and the self-similar 
lognormal solution of the corresponding Fokker-Planck equation. A very good agreement between 
the computed Boltzmann solution and the Fokker-Planck solution is observed. 



In the second test case we considered the most interesting situation with the presence of funda- 
mentalists, i.e both chartists and fundamentalists interact in the stock market. We compute an 
equilibrium situation where pf = Pc — 0.5 and the price stationary at the fundamental value 
Sf = 20. We take (3 — 0.1, tc = 1, 7 = 1-3, ol\ — ot2 = 0.01. We report the result of the simula- 
tion for the price distribution at the stationary state. In Figure [5] we show the price distribution 




Test 1 



Test 2 
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Figure 5: Stationary price distribution for the price with pp — pc — 0.5. Figure on the right is in 
log-log scale. The continuous line is the solution of the corresponding Fokker-Planck equation. 

together with the steady state of the corresponding Fokker-Planck equation. The emergence of a 
power law is clear also for the Boltzmann model, and deviations of the two models is observed for 
small values of the price. 

Test 3 

In the third test we consider the case with strategy exchange between the two populations of 
interacting agents. The switching rate used to run the simulation has the following form 

B FC (x) = e ax , 

where a represent the inertia of the reaction to profit differentials. We start the simulation con- 
sidering pc = Pf = 0.5. The fundamental price is Sp — 20, we take j3 = 6, tc = 0.02, 7 = 0.1, 
a = 0.8, fi = 0.2, D = 0.004, k = 0.75, and ip(y) — sgn(y). Furthermore we consider an herding 
function of the form H(y) = (1 — \y\). We run different simulations for T — 2000 iterations, with 
different values of a\, and 012, which measures respectively the herding and the market influence 
on the chartists. Three fundamental behaviors can be observed. The predominance of chartists, 
which leads the market towards a crash or a boom (see Figure 15]), the predominance of funda- 
mentalists, which originates damped oscillation of the price towards the fundamental value (see 
Figure [7]), and a balanced behavior, characterized by periods with oscillation of the price around 
the fundamental value (see Figures [5] and • From the simulations it is observed that, if we start 
with a balanced population between chartists and fundamentalists, the parameter «2, which char- 
acterize the influence of the price trend on the chartists investment propensity, plays a determinant 
role in the competition between the two different trading strategies. In particular when a,i > 0.6 
fundamentalists are predominant and price oscillations become dumped. 

5 Conclusion 

We derived an interacting agents kinetic model for a simple stock market characterized by two 
different market strategies, chartists and fundamentalists. The kinetic system couples a description 
for the propensity to invest of chartists and the price formation mechanism. The model is able to 
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Figure 6: Market crash due to a chartist predominance. The plot has been magnified to keep 
the price scale constant. The chartist dynamic is characterized by the parameters a\ = 0.2 and 
a2 = 0.55. Figure on the left represent the price averaged over N s — 50000 samples. Figure on 
the right represent the variation of the chartists 's fraction among the entire population of agents. 
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Figure 7: Dumped oscillation in the mean price, due to a predominance of fundamentalists. The 
chartist dynamic is characterized by the parameters a\ = 0.2 and ai = 0.7. Figure on the left 
represent the price averaged over N s = 50000 samples. Figure on the right represent the variation 
of the chartists's fraction among the entire population of agents. 

describe several market phenomena like the presence of of booms, crashes, and cyclic oscillations of 
the market price. The equilibrium behavior has been studied in a suitable asymptotic regime which 
originates a system of Fokker-Planck equation for the chartist's opinion dynamics and the price 
formation. We found that in a system of agents acting only using a chartist strategy the distribution 
of price converges towards a lognormal distribution. This is in good agreement with what previously 
found in [5] and observed in [TB]. When a second strategy based on a fundamentalist approach 
is introduced in the model the prices distribution displays Pareto power law tails, which is in 
accordance to what observed in the real market data. In the description of the chartists behavior 
we also introduced a value function which takes into account the effect of some psychological 
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Figure 8: Oscillations with different amplitudes in the mean price. The price is computed averaging 
over N s = 500 samples. The chartist dynamic is characterized by the parameters a\ = 0.5 and 
Q!2 = 0.4. Figure on the right represent the variation of the chartists's fraction among the entire 
population of agents. 




factors in the opinion formation dynamic. The main effect is to introduce market instabilities and 
to reduce the number of stable equilibrium configurations of the system. Let us finally conclude 
by observing that in principle several generalizations are possible. We mention here the possibility 
to include multiple interacting strategies and/or the influence of the wealth as an independent 
variable in the market dynamics. 

A Fokker-Planck asymptotics for the agents distribution 

We report in this appendix the details of the derivation of the Fokker-Planck equation (|15aj) for 
the distribution of chartists. Following [35] first we recall the definition of weak solution for kinetic 
equations of the form © and ©. Let / = [-1, 1] and M p {I) — {(9 € M p : fj \y\ p dO(y) < +00} 
be the space of all Borel measure of finite p-th order momentum, equipped with the topology of 
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weak convergence of the measures. Let F s (I) be the class of all real functions h on I such that 
/i(±l) = h'(±l) = and hS m \y) is Holder continuous of order 5 

ii ft <->ii,= 3 „p ""'■■'-';"■'' <„ ( 2 o) 

yi^V2 12/1 — 2/2 r 
where < <5 < 1, m + 8 = s and denotes the m-th derivative of h. 

Definition 1 Let f a (y) € M p (I) with p > 1 an initial probability density, a weak solution for f3J) 
is any probability density f € C 1 (K + , M p (I)) satisfying 

4 / f(y,t)<f>(y)dy = [ [ B(y,y*,r),r)*)f(y)f(y*)(((>(y') - <p(y))dr}dr)*dy*dy (21) 



dt 



/or £ > and all (f> € F p (I), and such that 

lim / f{y,t)<j>{y)dy = I f{y)cj>{y)dy. 



t->o > 

The scaled density f(y, t) defined in (fl2|) satisfies the equation in weak form 

^ jj{y,T)<t>{y)d y = IJ^ J Q(vMv*)f(y)f(y*my') - <t>(y))dvdv*dy,dy, (22) 

where J C R is a suitable symmetric support for the random variable 77 which avoids the dependence 
of the kernel B on the variables y and y* . 

Given 5 > let us take <f> € F 2+ s(I). 
From the microscopic dynamic of chartists we have 

y' -y= a 1 H(y)(y - y») + a 2 ($ - y) + D(y)r). 

In the asymptotic limit £ — > 0, <r 2 — > 0, we have y — y' ^ and we can use the Taylor expansion 

0(2/) - <Kv) = (aiH(y)(y - y*) + « 2 (* - y) + D{y)v) 4>'{y) 

- (aiif(j/)(y - y,) + a 2 (l> - y) + D{y)^ <j>"{y), 

where, for some < 8 < 1 

y = + (1 - %. 

Inserting this expansion in the weak formulation of the Boltzman equation, we get 
— I f{y,T)(/)(y)dy = 

e(ij)e(-),) [{a x H{y){y - y*) + a 2 (l> - y) + D(y)r?) 4>'{y) 

Hy)Ky*)dridri*dy*dy 



+ 2 



dr 
1 



- (a x H{y){y - y.) + a 2 ($ - y) + D(y)r)) <l>"(y) 



2 



where 



Rfov) = ^J i J j2 ®(v)®(v*)(aiH(y)(y-y*) + a 2 ($-y)+D(y)T i y 
■ 0"(y) - <f>"(y))f(y)f(y*)dri*dridy*dy. 



(23) 
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In order to prove that the remainder goes to zero as £ — > we start observing that, being 
<f) £ F 2+ s(T), and \y - y\ = 0\y' - y\ we get 



Hence 



\R(Z,*)\ < 



- 4>"(y)\ < \\4>"h\y - v\ 5 < WW - y\ 5 - 

Wh - - 

2£ 

2+5 _ 

+ D(y)r) 
pplying 1 

2+6 



6(77)6(77,) • 

aiH(y)(y - y*) + a 2 ($ - y) + D(y)-q f(y)f(y^)dn^dr]dy^dy. 
Using the fact that \H(y)\ < 1, |$| < 1, \y\ < 1 and applying the following simple inequality 

ai H(y)(y - yj + a 2 (£ - y) + £(77)77 ~ ' ' < C s {a 2 + 5 + a 2+s + \r,\ 2+s ) 
with Cs a suitable positive constant, we finally obtain 



\Rfo*)\ <c sP 2 c^ 



a 2+6 + a 2+s + I Q^M^dq 



\2+5, 



To simplify computations, we assume that 6, with zero mean and variance A£, is the density of 
i/X^W , where W is a random variable with zero mean and unit variance, that belongs to M 2 _|_ a , 
for a > 5, so we have 



J ®(rj)\v\ 2+S dr] = E 



2+8 



2+S\ 



and E (|W| 2+<5 ) is bounded. This is enough to show that in the asymptotic limit defined by (fT5| 
the quantity i?(£, a) tends to zero. 



Finally taking the limit in the weak formulation yields 
1 



lim 



i 2 J j 2 



6(77)6(77,) a x H{y){y - y.) + a 2 ($ - y) + D(y) V ) </>'(y) 



+ \ (aiH{y)(y-y*) + a2($-y)+D(y)r)) <f>"{y) 



f(y)f(y*)dr]dri*dy*dy 



( 



pcdxHWOT -y) + Pc a 2 ($ - y) ) 4>'{y) + '-( Pc D 2 {y))4>'- \y) 



A 



f{y)dy, 



which is nothing but the weak form of the Fokker-Planck equation (|15a[) . We can then state the 
following theorem 

Theorem 1 Let the probability density f £ Mq{I), and let the symmetric density 6 be in M 2+a 
with a > 8. Then in the asymptotic limit defined by H13\) the weak solution to the Boltzmann 
equation A22\) for the scaled density f{y, r) converges, up to extraction of a subsequence, to the 
weak solution of the Fokker-Planck equation \15a\) . 



B Fokker-Planck asymptotics for the price distribution 

In this appendix we derive the Fokker-Planck limit (| 15b[) for the scaled density distribution of the 
price. Now let .F S (R + ) be the class of all real functions ZionK such that h(0) — h'(0) = and 
h m (y) is Holder continuous of order S. We have the following 
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Definition 2 Given an initial price distribution Vq(s) G M p 
(0) is any probability density V £ C 1 (1R + , M p (M. + )) satisfying 



with p > 1 a weak solution to 



- I V(s,t)cf>(s)ds 



b(s,r))V{s,t)((j)(s') - (f)(s))dr)ds 



(24) 



for t > and all <f> G i 7 'p(R + ) and such that 

V(s,t)<p{s)ds 



lim 



4>(s)V (s)ds. 



Again we start with the weak formulation which now reads 



dr 



— I V(s,T)<t>( S )ds = ~ 



+ Jk 



^(r])V(s)(<f>(s') - (j)(s))dr,ds, 



(25) 



where K C K is a suitable symmetric support for the random variable r\ which avoids the depen- 
dence of the kernel b on the variable s. 

Let us take (j> £ i^+^R 4 ") with S > 0. Using a Taylor expansion of 4> around s 

0( s ')_0( s ) = (p(p c t c Ys + p F j(S F -s)) + r,s)<t>'(s) 



+ ~ (P(pct C Ys + p Fl {S F - a)) + ns) 2 <f>"(a), 



where for some < 9 < 1 



.5 = 9a' + (1 



and substituting into (|25[) we have 



— / V(s,T)</>( S )ds 



+ Jk 



^(vW(pctcYs + PFl {S F - a)) + ns) cj>'(s) 



where 



+ Jk 



+ 2 (PiPctcYs + PFl {S F - a)) + rjsY ct>"(a)]V(a)dr)ds 



(KpctcYs + PFl {S F - a)) + ns) 2 • - 4>" {s))V {s)dnds. 



Analogously as before, in order to perform the asymptotic limit we need to show that the 
quantity R((3, £, £ ) approaches zero as £ — >■ 0. We observe that being g F 2+( s(IR+) and |s — s| = 



s\ we have 



\ct>"(S)-cj ) "(s)\ < \\d>"\\ s \ S ' 



hence 



w,c,o\ < 



R+ JK 



PlpctcY 



Pf7- 



(S F -s) 



2 + ,S 



„2+S 



V(s)drjds 



Next we observe that 



P (pctcY + Pf1 {Sf s S ' 
C 2+S ({!3pctc) 2+s + {Ppfi) 2+s 



2 + ,S 



< 



S 2+S + S 2+S 
Q 2+5 



(26) 



\2+5 
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where C2+8 > is a suitable constant. 

As in appendix A we assume that ty, with zero mean and variance is the density of \fvCyV , 
where W is a random variable with zero mean and unit variance, that belongs to M2+ Q , for a > 8, 
so we have 

'uCW 



J ^(v)M 2+S dv = E (\ \Z^W\ 2+S ^j = (vC) 1+ ^E (\W\ 2+S ) , 



(27) 



and E(\W\ 2+S ) is bounded. 
Then we obtain 



W.Cfll < C 2+5 Mi{[(/3^c) 2+5 + (^7) 2+5 + K) 1+ ^(l^| 2+5 ) 



2i 



„2+5 



y( s )d s + (f3 PFl ) 2+s s 2 F +s 



From this inequality it follows that R((3, £, £) tends to zero in the limit (| 14[) if 

l/(s,r)s 2+5 ds 



is bounded at any fixed time r > 0, provided that the same bound holds at time r = . 
To show this we start again from the weak formulation ([25]) . The choice <p(y) — y p gives 



d_ 

dr 



V(s)s p ds = - 



R+ J K 



V(r])V(s){s' p - s p )dr)ds. 



Now 



s'p - s p= psP- 1 ^' -s) + -p(p - l)s p - 2 (s' - sf 



where for some < 6 < 1 

s = 0s' + (1 - 6)s. 

Recalling the microscopic dynamic for the evolution of the price variable s we can write 



d_ 



V(s)s p ds 



1 

£ Jm+ jk 
P 



p S P-\s' - s) + -p{p - 1)~ S p- 2 {s' - s) 2 



drjds 



+ Jk 



2£ 



^(t 1 )V(s) S p- 1 [(/3{p c t c Ys + p Fl {S F - s)) + V s} dr,ds 



V(ri)V(s)s p - 2 [(3 (pct C Ys + pfi{S f - s)) + tjs} 2 drjds 
+ Jk 



Since the random variable rj has zero mean value, the first term in the last expression reduces to 



(pctc'Y-PFj) V(s)s p ds + ppSpl V(s)s p - 1 ds 



For the second therm, we know that 



s + l3(p c t c Ys + p Fl {S F - s)) + rjs) + (1 

(s F - s y 



PctcY + PFl- 



6r)+l 
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which implies 



qP- 
o p 



-2 + gP -2 



sP- 



+ i?7r 2 + i 



with C p a suitable constant. 



Gathering all this the weak formulation gives 





d_ 



V(s)s p ds < p- 



(pctcY-ppl) / V(s)s p ds + pfSfI / V(s)s p - 1 ds 



P(P ~ !) n 



B+ J K 



^{n)V{s)s p 



P PctcY + PFl 



(S F - s) 



{ppctcY' 2 + {PpfiY' 2 



q p-2 



S p ~ 



IP-2 , 1 



drjds. 



Now if we consider the asymptotic limit (jl4j) and recall (|27|) for the high order moments of n, it 
follows that the p-moments of V(s,t) are bounded at any finite time independently of £ and for 
p > 2 + 5 satisfy 

— / V{s)s p ds<A p [ V(s,t)s p ds + B p [ V(s,t)s p ~ 1 ds 
dr J ffi + J R + J R + 

where A p — p(3 (pctcY — PFl) + p(p — t)vC p /2 and B p = p(3p F SFl- 
Coming back to the asymptotic expansion we can finally perform the limit 



lim — 



R+ JK 
+ 



*(V) W(p c (t)Yt c s + P Fl(S F - s)) + vs) 0(a) 

V(s)dr/ds 



i (p{p c (t)Yt c s + P fi(Sf - a)) + f]s) 2 cf>"(s) 



P(pc(t)Yt cS pFl(S F ~ s))<t>'(s) + -sV(s)J V(s)ds, 

which is the weak form of the Fokker-Planck equation for the price (|15b[) . So we proved the 
following 



Theorem 2 Let the probability density Vq G Mo(R + ). Then in the limit defined by {1J$ the weak 
solution to the Boltzmann equation i25\) for the scaled density V(s, r) converges, up to extraction 
of a subsequence, to a weak solution of (| 15b\ . 
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